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We revisit sneutrino dark matter in light of the recent results from the PAMELA and ATIC 
- . . experiments. In the U{1)b-l extension of the minimal supersymmetric standard model the right- 

0^ , handed sneutrino is a natural candidate for thermal dark matter. Sneutrino annihilation at the 

^— . present time can be considerably enhanced due to the exchange of the lightest field in the Higgs 

»' ' sector that breaks U{1)b-l- The annihilation mainly produces taus (or muons) by the virtue of 

B — L charges. A sneutrino mass in the 1 — 2 TeV range provides a good fit to the PAMELA 
O . data and a reasonable fit to the ATIC data. Within this mass range the sneutrino-nucleon elastic 

D ' scattering cross section is 10~^^ — 10~^ pb, which might be probed by upcoming and future direct 

pL^ I detection experiments. In addition, if (at least) one of the neutrinos is dominantly a Dirac fermion, 

, the sneutrino can provide a unified picture of dark matter and inflation. 
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I. INTRODUCTION 

Even though the existence of dark matter has been sup- 
ported by various lines of evidence, the identity of dark 
matter itself is not yet known. One proposed solution 
for this dark matter problem comes from particle physics 
beyond the standard model in the form of weakly inter- 
acting massive particles (WIMPs) In particular, for 
weak scale masses and interactions, thermal freeze out 
of WIMP annihilation in the early universe can result 
in an acceptable relic abundance for dark matter, as pre- 
cisely measured by cosmic microwave background (CMB) 
experiments Supersymmetry, as one candidate for 
physics beyond the standard model, has a natural dark 
matter candidate in the lightest supersymmetric particle 
(LSP). It is known that in supersymmetric models a neu- 
tralino LSP is a suitable candidate for dark matter [3]. 

There are currently major experimental efforts for both 
direct and indirect detection of the dark matter particle. 
Direct detection probes the scattering of the dark matter 
particle off nuclei in underground dark matter detectors, 
while indirect detection investigates astrophysical effects 
of dark matter annihilation in the galaxy, including signa- 
tures in the cosmic rays. PAMELA is a satellite-borne ex- 
periment that measures cosmic ray fluxes. The recently 
published results show an excess of positron flux at en- 
ergies above 10 GeV (J], while no excess of anti-proton 
flux is observed The publication shows results up to 
~ 100 GeV and the experiment is expected to get data 
up to ^ 190 GeV for anti-protons and ~ 270 GeV for 
positrons. Another cosmic ray experiment called ATIC 
(a balloon experiment) has also recently published data 
where one observes an excess in the e"*" -I- e~ spectrum 
with a peak around 600 GeV . There is a third exper- 
iment, the PPB-BETS Q (also a balloon experiment), 
which reports an excess in the + e~ energy spectrum 
between 500 and 800 GeV. However, the excess is based 
on a few data-points that are not quite consistent with 
the ATIC data. While there could be astrophysical expla- 



nations for these anomalies (e.g. from nearby pulsars [Sj]), 
it is reasonable to ask whether they can be attributed to 
the effect of dark matter annihilation in the galaxy. 

Model-independent analysis shows that the annihila- 
tion cross section required to explain the positron ex- 
cess exceeds the canonical value required by relic density, 
i.e. '--^ 3 X 10~^^ cm^/s, by at least an order of magni- 
tude 0. In the usual neutralino dark matter scenario 
in the minimal Supergravity (mSUGRA) model, the sit- 
uation is further complicated because the dark matter 
annihilation (to fermions) in that model is P-wave sup- 
pressed, implying a much smaller annihilation cross sec- 
tion today as compared to that at the freeze out time. 
An astrophysical boost factor of lO'^ — 10'' js then needed 
to explain the observed positron excess [iq|. However, 
this might be difficult to obtain based on recent anal- 
yses of halo substructures (see e.g. [HI). Moreover, in 
order to explain both the positron and anti-proton data, 
dark matter annihilation must be dominated by leptonic 
final state modes [12, [3 • (There could also be some ef- 
fects from anisotropic propagation on the positron and 
anti-proton fluxes that still need to be investigated [l^.) 
There have been proposals [HjUBI (also see [Ifi^) for new 
dark matter models in which the dark matter candidate 
belongs to a hidden sector, and an acceptable thermal 
relic density is obtained via new gauge interactions. The 
key ideas of these models are that the dark matter an- 
nihilation today is enhanced by a Sommerfeld effect [l3| 
due to the existence of light bosons and that annihila- 
tion mainly produces lepton final states via symmetry of 
the hidden sector. This arrangement explains PAMELA 
data for a dark matter mass of a few hundred GeV, with- 
out needing a large astrophysical boost factor, and ATIC 
data for larger values of dark matter mass 'la. Another 
type of explanation that has been proposed for the data 
is decaying dark matter with a tuned lifetime p^ . 

We recently proposed an explicit model that can ex- 
plain the measured anomalies in the cosmic rays [20| . It 
is based on a simple extension of the minimal supersym- 
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metric standard model (MSSM) that includes a gauged 
U{1)b-l and where the dark matter is the lightest neu- 
tralino in the new sector. Even though this model has a 
large dark matter annihilation cross section today due to 
Sommerfeld enhancement, the cross section for scatter- 
ing of dark matter off quarks is too low to be accessible 
to direct detection experiments. In fact, this is a generic 
situation for hidden sector dark matter models that can 
explain PAMELA data along the line discussed above. 

In this paper we again consider a B — L extension of 
the MSSM, but with the right-handed (RH) sneutrino 
as the dark matter. As we will argue, this is a mini- 
mal model of thermal dark matter that can explain the 
observed anomalies in the cosmic rays and can also be 
probed by direct detection experiments. The main chan- 
nel of sneutrino annihilation is to light Higgs fields, which 
carry a non-zero B -- L quantum number. These Higgs 
particles in turn decay dominantly to leptons by virtue of 
the B — L charges for fermions. The same Higgs field also 
results in a large Sommerfeld enhancement factor for the 
annihilation cross section. For a sneutrino mass of 1-2 
TeV, this model can explain PAMELA and ATIC data. 
In addition, due to the scalar nature of dark matter, the 
sneutrino-proton elastic scattering cross section is in the 
10~^^ — 10~^ pb range, which is an interesting range from 
a direct detection perspective. Moreover, the sneutrino 
can be part of the field that drives primordial inflation, 
thus explaining the small temperature anisotropy in the 
cosmic microwave background (CMB) via tiny neutrino 
masses [2l|, . We will also discuss various possibilities 
for radiative breaking of the B ~ L symmetry and some 
related issues. 



II. THE MODEL 

The B ~ L extension of the MSSM [H] is well moti- 
vated since it automatically implies the existence of three 
RH neutrinos through which one can explain the neu- 
trino masses and mixings. The minimal model contains 
a new gauge boson Z', two new Higgs fields H[ and H2, 
the RH neutrinos N, and their supersymmetric partners. 
The superpotential is (the boldface characters denote su- 
perfields) 



MSSM 



Wb-l + , 



(1) 



where Hu and L are the superfields containing the Higgs 
field that gives mass to up-type quarks and the left- 
handed (LH) leptons respectively. For simplicity, we have 
omitted the family indices. The Wb~l term contains 
H']^, H2 and N'^. Its detailed form depends on the charge 
assignments of the new Higgs fields (explicit examples 
will be presented later). The last term on the RH side of 
Eq. ID) is the neutrino Yukawa coupling term. 

The scalar potential consists of F-terms from the su- 
perpotential, and D-terms from the gauge symmetries. 



The Z^-term contribution from U{1)b-l is given by 

1 



where 



Db-l — 29b-l 



Qii\H[f-m') + l\N\' 



(2) 



(3) 



Here gs^L is the gauge couphng of U{1)b-l, and +Qi, 
—Qi, 1/2 are the -B—L charges of H2, respectively 
{N is the sneutrino field). The U{1)b-l is broken by the 
vacuum expectation value (VEV) of H[ and H2, which 
we denote by v[ and V2 respectively. This results in a 
mass mz' — gB-hQix/v'i + v'2 for the Z' gauge boson. 
We have three physical Higgs fields (/>, $ (scalars) and A 
(a pseudo scalar). The scalar Higgses are related to the 
real parts of H2 through the mixing angle a': 



v'l + cos a'<^ — sin a' (j) H'u 



V2 



V2 



I v'2 + sin + cos ii'2 I 
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V2 



V2 ' 



(4) 



where H[j,H2j represent the imaginary parts. 
Eqs. (|2l3l4p lead to the following terms in the scalar 
potential 

-^.9L.QiCOs(2a')0'|iV|' 

1 ~ 2 

+ ^gB-Lmz' cos(a' + (3') $ \N\ 

+ ... , (5) 

where tan/3' = v'2/v'i. The masses of the Higgs fields 
follow < cos^(2/3') m\, and m$, ^ mz>. 

A natural dark matter candidate in this model is the 
sneutrino TV ^. We note that the N has fewer gauge 
interactions than other fields, hence its mass receives 
the smallest contribution from the gaugino loops. The 
main processes for annihilation of dark matter quanta 
are then governed by interactions in Eq. ([5]). The dom- 
inant channel is N*N cfxj) via the s-channel exchange 
of the (p, the t, u-channel exchange of the A'', and 
the contact term \N\'^(j)^. The s-channel Z' exchange is 
subdominant because of the large Z' mass (as required 
by the experimental bound on mz'). There are also 
N*N 0$, (f>A, <f><I', AA annihilation processes, but 
they are kinematically suppressed and/or forbidden for 
the parameter space we are considering. The sneutrinos 



^ Another candidate is the hghtest neutralino in the new sector, 
which is a linear combination of the U{1)b—l gaugino Z' and 
the two Higgsinos H[. Hi, [20l[2^. 
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can also annihilate to RH neutrinos via f-channel neu- 
tralino exchange. Again for the parameter space that we 
consider the annihilation into vD final states is at least 
an order of magnitude below the 00 final states. Other 
fermion final states, through s-channel Z' exchange, have 
even smaller branching ratios. Moreover, note that the 
annihilations to fcrmion-antifcrmion final states are P- 
wave suppressed. 

The (j) subsequently decays into fermion-antifermion 
pairs via a one-loop diagram containing two Z' bosons. 
The decay rate is given by: 

(6) 

where Q / and Q^p are the B — L charges of the final state 
fermion and the respectively, is the fermion mass, 
and Cf denotes color factor. Since the B — L charge 
of leptons is three times larger than that of quarks, the 
leptonic branching ratio is naturally larger than that for 
quarks. We note that can be controlled by the VEVs 
of the new Higgs fields and for comparable VEVs, i.e. for 
tan /J' « 1, it can be very small compared to mz'. For 
> 2mb the dominant decay mode is — *■ t^t+ final 
state, while the branching ratio for the (j) ^ bb mode is 
« 7 times smaller. 

The annihilation cross section at the present time has 
Sommerfeld enhancement as a result of the attractive 
force between sneutrinos due to the exchange. The 
Higgs coupling to dark matter is given by the first term on 
the RH side of Eq. ^ and leads to an attractive potential 
V{r) — —a{e~^*^/r) in the non-relativistic limit [13], 
where 

gB-L-mz' sin(a' + /?') 
a = z : (7) 

and m>7 is the sneutrino mass. 



III. SNEUTRINO DARK MATTER AND 
PAMELA 

We are now going to show that the sneutrino dark mat- 
ter can explain the PAMELA data. We first identify the 
allowed regions of the model parameter space that re- 
sult in an acceptable dark matter relic density and then 
find the Sommerfeld enhancement factor for these re- 
gions. For an explicit example, which we call the minimal 
model, we choose the B — L charge for H[ (i.e. Qi) to be 
3/2. The B — L charges of the fields involved are shown in 
Table HI We use reasonable values for the model param- 
eters, i.e., tan/3' « 1, mz' > 1.5 TeV, ^i' = 0.5-1.5 TeV 
{fi' being the Higgs mixing parameter in the B—L sector), 
soft masses for the Higgs fields niH' , = 200 — 600 GeV, 
and soft gaugino mass M^, > 1 TeV. We use gg-L ^ 0.40, 
which is in concordance with unification of the gauge cou- 
plings (we need to use a normalization factor \/3/2 for 



Fields 


Q 




L 




H[ 




Qb-l 


1/6 


-1/6 


-1/2 


1/2 


3/2 


-3/2 



TABLE I: The B — L charges of the fields for the minimal 
model. Here Q and L represent quarks and leptons respec- 
tively, while H[ and H2 are the two new Higgs fields. The 
MSSM Higgs fields have zero B — L charges. 

unification). We show the unification of all the gauge 
couplings using the two loop renormalization group equa- 
tions (RGEs) in Figure [1] We find that for m^/ ~ 2.5 
TeV the couplings unify at ~ 10^^ GeV. This figure is 
drawn for the B — L assignments shown in Table [J The 
Z' mass used in the calculation obeys the LEP and Teva- 
tron bounds [1^ |2^ for our charge assignments. The 
sneutrino mass is chosen to be between 800 GeV and 2 
TeV in order to explain the PAMELA (and ATIC) data. 
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FIG. 1: We show the unification of gauge couplings for the 
B — L charge assignments in Table 1 using two loop RGEs. 

In Figure [H we show the relic density and the sneu- 
trino mass for different model points. The horizontal 
band shows the acceptable range for relic density accord- 
ing to the latest WMAP data ^] . In Figure[3]we show the 
possible Sommerfeld enhancement factor R that can be 
obtained for these points in term of = m^/amfj. Note 
that many of the points that satisfy the relic density con- 
straint have enhancement factor R > 10'^, corresponding 
to = 0.55 to 0.65. This is true for the whole range of 
TOjY shown in Figure O The lifetime of 4> for these points 
is found to be ^ 10~^ — 10~^ seconds from Eq. ([6]). 
Thus (p's produced in the early universe decay rapidly 
enough in order not to affect big bang nucleosynthesis 
(BBN) [27] 2. We select points that satisfy the dark 
matter relic density and then use DarkSUSY-5.0.2 [2§] to 



^ Since dark matter particles are non-relativistic at the time of 
BBN, their annihilation enhanced by Sommerfeld effect can re- 
sult in significant electromagnetic and/or hadronic showers that 
dissociate light elements from BBN. For muon final states, the 
large enhancement factor required to explain PAMELA is com- 
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FIG. 2: We show the rehc density and the sneutrino mass for 
model points generated by varying the parameters mentioned 
in the text. 




I 



4000 

3000 

2000 

1000 
500 
100 



FIG. 3: We show rehc density as a function of e^. We show 
different ranges for the Sommerfeld enhancement factor R by 
shaded contours. 



up to ~ 15 GeV and still have acceptable anti-proton 
flux. We can also have a reasonable fit to the ATIC data, 
although simultaneous fit for both ATIC and PAMELA 
are not satisfactory j20j. 

We note that this model also has a great potential to 
be observed with the Fermi Satellite experiment. Due 
to electromagnetically charged final states of (p decays, 
the Sommerfeld enhancement would also lead to a higher 
rate of photons in the gamma ray background [33j|. There 
could also be contribution to the neutrino flux 1341. 
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FIG. 4: We show a fit to the PAMELA data when the </> 
decays mostly to taus (dark blue) or muons (light green) for 
a sneutrino mass of 1.5 TeV and an enhancement factor of 
10^. The dashed line is the expected background cosmic rays. 



calculate the positron flux from dark matter annihila- 
tion. Each pair annihilation in our model produces 2 0's 
that yield four fermions upon their decay. For this rea- 
son, we generally need a heavier sneutrino compared to 
models in which the pair annihilation directly produces 
fermions. We normalize the positron fraction by a fac- 
tor kb — 1.11 [30]. There are theoretical uncertainties 
in the positron cosmic ray flux calculation due to the 
assumptions about the dark matter halo profile and the 
cosmic ray propagation model. Here we assume NFW 
profile [3l| for the dark matter halo and MED parame- 
ters for the propagation as defined in [33 |. 

In Figure m we show our fit to the PAMELA data for 
= 1.5 TeV for r+r" and ijl^ final state cases. We 
found that with an enhancement factor of lO'^ the chi- 
square values (including only points with energy greater 
than 10 GeV) for a sneutrino mass of 1.5 TeV are small, 
i.e. 2.9 and 5.5 for T+r~ and jX^ IJ.~ respectively. When 
is (chosen to be) below 2mf, but above 2mr, we do 
not have any anti-proton excess. In fact we can raise rriA, 



patible with BBN bounds, while for tau final states a small as- 
trophysical boost factor seems to be needed in order not to affect 

BBN m. 



IV. DIRECT DETECTION 

The current upper bound on the spin-independent 
dark matter particle-proton scattering cross section is 
about 4.6 X 10^* pb for a dark matter mass around 
60 GeV, and increasing to ~ 2 x 10"^ pb for a mass 
around 1.2 TeV [35] ■ In our model the elastic scatter- 
ing of the sneutrino occurs via the Z' exchange with 
the nucleus in the t-channel. This leads to only a spin- 
independent contribution since the B — L charges of 
the left and right quarks are the same. In Figure [5l 
we show the N-p scattering cross section for the model 
points in Figure 1 that satisfy the relic density constraint 
0.096 < ^DAih^ < 0.124. We see that the cross section 
can be in the 10^^^ — 10~^ pb range, which is close to 
the reach of the upcoming dark matter direct detection 
experiments [36j . 

It is also seen that the cross section decreases as the 
sneutrino mass increases. This is because for larger val- 
ues of TO^ we also need a larger annihilation cross section 
to satisfy the relic density constraint. As discussed earlier 
the annihilation cross section depends on the sneutrino 
couplings to (j) and see Eq. ([5]), which are oc mz'- 

It is interesting to note that within this mass range 
the Z' can be produced at the LHC. The Z' decay will 
produce the new light Higgs (among other fields). How- 
ever, (j> will decay outside of the detector because of its 
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FIG. 5: We show the direct detection cross section as a func- 
tion of sneutrino mass. 



relatively long life time (~ 10^^ sec). Thus, in addition 
to the sneutrino LSP, we have another source of missing 
energy signal in this model ^. 



V. RADIATIVE BREAKING OF B L 
SYMMETRY 

In order to have spontaneous breaking of U{1)b-l we 
need a negative eigenvalue in the H[, H2 square-mass 
matrix. This can arise dynamically as a result of radia- 
tive corrections to the Higgs soft masses. Here we discuss 
two generic possibilities. 



A. Higgs coupling to right-handed neutrinos 

If Qi = 1, then H2 can have a superpotential coupling 
to the RH neutrinos In this case we have: 



(8) 



where /z' is analogous to the MSSM /i parameter. Note 
that in the minimal model, with Qi — 3/2, the H2 can 
only couple to H[ in the superpotential. Taking the soft 
mass parameters into account, the Higgs potential is 



{Bn'H[H'2 + h.c) 



(9) 

(Here the parameters mi,m2,B are not to be confused 
with the MSSM Higgs parameters.) The Yukawa cou- 
pling between and N'^ can drive m| to a sufficiently 



negative value such that + n'^ < around the TeV en- 
ergy scale ^. This requires that / not be very small. On 
the other hand, / should not be very large. Otherwise 
there would be a one-loop correction that lifts the 4> mass 
above its tree-level limit m| < Af|,cos^(2/3') (similar to 
the correction from the top Yukawa in the MSSM [ssj): 



At 



(10) 



Here tun — is the Majorana mass of N and 



"^Nr = [m\ + mN{A + ^')\ 



1/2 



1/2 



(11) 



are the total masses of the real and imaginary parts of 
the N respectively. The mass splitting is due to the con- 
tribution of the A and ji' terms to the N potential given 
by 



f{AH'2 + fi' H[*) NN + h.c. 



(12) 



For simplicity we have assumed that /, A and /i' are 
real, and we have used the requirement that v'l ~ V2. 
This mass splitting can lift above its tree-level limit 
similar to the effect from stop mixing in the MSSM [38| . 
In order for this not to be significant, it will be sufficient 
to have |A -I- <C m^y at the TeV scale ^ . 

In our case the sneutrino and the neutrino masses can 
be close and therefore the corrections to the Higgs mass 
can be very small < 10 GeV (which we have checked nu- 
merically). In fact, this allows us to keep the Higgs mass 
below 15 GeV such that we do not have any problems 
with anti-proton data. We also have found that / ~ 0.2 
is large enough to lead to radiative breaking of U (1)_b-l- 
Assuming grand unification, one should use the SUGRA 
boundary conditions to achieve radiative breaking due 
to the smallness of / (as compared with the top Yukawa 
in the MSSM) . The appropriate boundary conditions are 



such that 



< 0, but 



n,2 + 



> 0, while all 



other soft (mass)^ are positive at the grand unification 
scale. 

The potential term in Eq. (|12p. no matter how small, 
inevitably splits the masses of the real and imaginary 
parts of the sneutrino. In this case dark matter will be 
the lighter of and iV/, which we choose to be Nr 
without loss of generality. Then the heavier component 
Nj decays to dark matter and fermions via an off-shell Z' 
in the early universe. In this case the dark matter cou- 
pling to the light Higgs (j), the first term on the RH side 



^ We also note that there are 7 neutralinos in this model, compared 
to four in the MSSM, while the number of charginos is still two. 
Therefore, using the end point analysis "st], one can find many 
neutral states. 

Similarly, H'^ can couple to the RH neutrinos if Qi = — 1. 



^ Since H'^ has no Yukawa couplings, increases towards smaller 

scales because of the U{1)b—l gaugino loop. 
^ This also ensures that the contribution of Eq. I I12II to the coupling 

of dark matter to <j>, <I> is negligible compared with that from the 

D-term JSJ. 
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of Eq. ^ , will be twice as small as the term wherein the 
dark matter is a complex field. However the attractive 
potential, due to (f) exchange, will be the same in the two 
cases. This is because a real field annihilates and creates 
the same quanta, thus giving rise to a factor of two that 
compensates the 1/2 that appears in the coupling. The 
same also holds for dark matter annihilation. Therefore, 
as far as relic density calculations and the Sommerfeld 
enhancement factor are concerned, our results will not 
depend on whether dark matter is the complex N field 
or its real (or imaginary) component. 

We also note that the /HaN'^N'^ term opens a new 
channel for the decay of the light Higgs 4>. Since <C 
TOjv, the decay to on-shell RH neutrinos is not possible. 
The decay can happen via off-shell N through their cou- 
plings to the MSSM fields ([1]), at higher orders of per- 
turbation theory. Then, if the neutrino Dirac Yukawa 
coupling is sufficiently small, this decay mode will be 
totally negligible. Indeed this is the case since the see- 
saw mechanism at the TeV scale requires yo <, 10~^ to 
generate the observed neutrino masses. 



B. Higgs coupling to other fields 

If Qi 7^ ±1, then neither H[ nor H2 can have a su- 
perpotential coupling to N. This happens, for example, 
if Qi = 1/2 [23] or Qi =3/2 (as chosen above for the 
minimal model and also in [20J). Radiative breaking of 
U{1)b-l then requires a moderate Yukawa coupling of 
one of the B — L Higgses to the new fields. 

A simple model of this type includes two new super- 
fields ^"1, 4>2 and has the following superpotential 



Wb-l = //h;h'2 

+ /i"*i*2 



A*lN'=N^ 



(13) 



where Qi = 2 and the B — L charges of $1, $2 are 
— 1, -1-1 respectively. All we need for radiative breaking 
of U{1)b^l is /i ^ 0.2, or /2 ~ 0.2 (one of the couplings 
can be very small) and SUGRA boundary conditions. 
This ensures sufficiently large loop corrections that drive 
the (mass)^ of the corresponding Higgs field to negative 
values around a TeV. In this case, since the Higgs fields 
are not coupled to N"^, there are no Majorana masses for 
the RH neutrinos, and hence neutrinos are of a Dirac na- 
ture. Also, there will be no splitting between the masses 
of the real and imaginary parts of the N, and thus the 
dark matter is a complex scalar field. 

Interactions in the first line of Eq. ((T3)) result in a 
2x2 mass matrix for $1, $2- For large enough fj," 
(i.e. /i" ~ v'l, V2) the mass eigenvalues will be larger 
than the sneutrino mass. Hence the corresponding mass 
eigenstates would quickly decay to the RH neutrino and 
sneutrino through the A^iN'^N'^ superpotential term 
in the early universe. The /iH'i4>i*i and A^iN'^N'^ 
terms also open a new decay channel for the light Higgs: 



(j> — > 4N via off-shell fermionic components of <I>i ^. How- 
ever, A can be chosen to be sufficiently small such that 
— > T~T"*" remains the dominant decay mode of 0. 

If Qi = 3/2, as chosen in Table 1, one can introduce 
four new superfields $1, $2, ^4 (with respective 
B-L charges -1, -1-1, -1/2, -I- 1/2) and the following 
superpotential 

Wb-l = //H;H^ + /iH;*i*3 + /2H^*2*4 
+ /i"*i*2 + m"'*3*4 + A*lN=N^ 

(14) 

Again the mass eigenvalues of the 4x4 mass matrix for 
$1, $2, ^3, ^4 can be made large enough such that the 
corresponding eigenstates rapidly decay to lighter fields 
via the A^iN'^N'^ superpotential terms. 

It is interesting to note that although we have used 
Qi =3/2 to obtain the results in the previous sections, 
they are largely independent from the exact charge as- 
signments of H[, This is because the major con- 
tributions to relic density calculations involve Z' mass in 
the coupling and Qi is already absorbed in the mass defi- 
nition, see Eq. ([5|). The direct detection cross section also 
remains unchanged for the same reason. The gauge cou- 
pling unification still occurs but requires a larger value 
of gB-L for IQil < 3/2. Therefore Eqs. (|8ll3ll4p can all 
yield thermal sneutrino dark matter with a large Som- 
merfeld enhancement factor and radiative breaking of 
U{1)b-l- 



VI. RIGHT-HANDED SNEUTRINO AND 
INFLATION 

In addition to being the dark matter candidate, the 
RH sneutrino can also drive inflation in the context of 
the U{1)b-l model [21], [22|. The gauge-invariant combi- 
nation N'^HuL forms a I?-flat direction under the whole 
gauge symmetry SU{3)c x SU{2)l x U{1)y x U{1)b-l- 
The flat direction field cp is defined as 



V3 ■ 



(15) 



The potential along the flat direction, after the minimiza- 
tion along the angular direction, is found to be [2l| . 



(16) 



where yjj is the neutrino Dirac Yukawa ([T|) and A is the 
corresponding j4-term coupling. The flat direction mass 
is given in terms of the N , L masses: 



m 



(17) 



^ Note that there can be no in the final state since <^ rrij-^. 



7 



For A w 4to;p, there exists an inflection point (po ioi 
which V"{(po) — 0. Due to the extreme flatness of the 
potential around the inflection point, inflation can take 
place near (po- The amplitude of density perturbations 
generated during inflation follows [2l[ 

(Im^)' (!;)^-.. (IS) 

where ~ 174 GeV and A/cobe ~ 50 is the num- 

ber of e-foldings between the time that relevant pertur- 
bations were produced and the end of inflation. It is 
seen from Eq. that perturbations of the correct size 
Sh = 1-91 X 10~^ are obtained if yo 10~^^. Interest- 
ingly this is the typical neutrino Dirac Yukawa coupling 
that gives rise to the mass scale ~ 0.05 eV required 
to explain the atmospheric neutrino oscillations detected 
by the Super-Kamiokande experiment [40| if neutrinos 
are dominantly Dirac in nature. 

Dirac neutrinos can be achieved in both cases consid- 
ered above (|8ll3p . If the Higgs fields that break U{1)b-l 
are not coupled to (|13p . the neutrinos are naturally 
Dirac since there will be no Majorana masses for the 
RH neutrinos in this case, regardless of how big the cou- 
plings /i, /2 are. On the other hand, a superpotential 
coupling between the Higgs that breaks U{1)b-l and 
N"^, see Eq. ([5]), inevitably induces a Majorana mass 
iriN = 2/(i?2) to the RH neutrinos upon the B — L break- 
ing. Nevertheless, the main contribution to the mass of 
light neutrinos comes from the Dirac Yukawa coupling 
Ud- Hence all that we need in this case is one of the Ma- 
jorana masses (out of the three generations) to be very 
small in order to have an almost Dirac neutrino with 
Ud ^ 10~^^. The other Majorana masses (and respec- 
tively the coupling /) can be large. 

We therefore conclude that it is possible to have a uni- 
fied U{1)b-l model of inflation and dark matter, where 
the RH sneutrino is the dark matter and a component of 
the inflaton field 1231. 



VII. CONCLUSION 

Motivated by the recently reported cosmic ray anoma- 
lies, we have reconsidered a minimal extension of the 
MSSM that includes a gauged U{1)b-l- This additional 
symmetry is broken around a TeV by two new Higgs 
fields that carry non-zero B — L charges. The RH sneu- 
trino can naturally be the dark matter candidate in this 



model since it has the smallest gauge interactions among 
all the fields. Sneutrino interactions of gauge strength 
yield an acceptable thermal relic density in large regions 
of the parameter space. If the lightest Higgs in the B — L 
sector is much lighter than a TeV, the dark matter dom- 
inantly annihilates into final states including this Higgs. 
The annihilation is governed by D-term couplings be- 
tween the sneutrino and the Higgs and takes place in the 
5-wave. The light Higgs subsequently decays to fermions 
and the B — L symmetry guarantees that the branching 
ratio for producing leptons is several times larger than 
that for quarks, which agrees with the observation of 
positron and anti-proton fluxes by PAMELA. For a 1 — 2 
TeV sneutrino, a Higgs mass < 15 GeV will result in a 
large Sommerfeld enhancement factor 0(10'^) in the an- 
nihilation cross section at the present time. This provides 
a good fit to the PAMELA data and a reasonable fit to 
the ATIC data. 

The sneutrino interacts with quarks via t-channel ex- 
change of the U{1)b-l gauge boson Z'. The interaction 
only has a spin-independent part since B — L is a vector 
symmetry. The resulting sneutrino-nucleon elastic scat- 
tering cross section is found to be 10""'^^ — 10~^ pb, which 
might be within the reach of future direct detection ex- 
periments. 

We have also discussed radiative breaking oi U{1)b~l- 
This requires that (one of) the B — L Higgs fields have a 
relatively large Yukawa coupling to the RH sneutrino or 
some other field. A Yukawa coupling ~ 0.2 is sufficient 
to induce radiative breaking while keeping the mass of 
the light Higgs low enough to give rise to considerable 
Sommerfeld enhancement of dark matter annihilation. 

Finally, if (at least) one of the neutrinos is dominantly 
a Dirac fermion, the sneutrino can be part of the infla- 
ton field in addition to being the dark matter. This is a 
very appealing scenario since direct and indirect detec- 
tion experiments not only probe dark matter in this case, 
but they can also reveal the interactions of the inflaton, 
which is supposed to be the most elusive particle in the 
universe. 
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